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1. Introduction 

> : 

In a seminal work about two decades ago, Witten proposed that strange quark matter 
(SQM), the combination of roughly equal number of up, down and strange quarks, might 
be the true ground state of quantum chromodynamics (QCD) [lj. Later calculations 
have shown that in spite of the effect of their finite volume, small nuggets of SQM in 
form of "strangelets" can also be stable as long as the baryon number exceeds a critical 
value A crit [2J. On the one hand, various theoretical scenarios have provided chances 
for strangelet formation [31 Hj. They could be produced in highly energetic nuclear 
collisions [5], might originate from the collisions of two strange stars [Sj, and could 
also be ejected by supernova explosions [7j. On the other hand, several exotic cosmic 
ray events have been reported by balloon and mountain experiments [8], [9], [10], which 
are considered to be ideal candidates of strangelets. The ultra-high energy cosmic ray 
(with energy > 10 19 eV ) events could also be the results of extensive air showers of 
relativistic strangelets accelerated in pulsar magnetospheres [3J, [TT| . Interestingly, one 
doubly charged event, with charge- mass ratio of ~ 0.1, has been detected by the AMS 
experiment in space [12] and suggested to be strangelet-originated. This idea could be 
tested in the future AMS02. 
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Since the existence of stable SQM would have remarkable consequence in 
cosmology and astrophysics, what is experimentally important is to find out strangelets' 
contribution to cosmic ray flux and the mechanism for the propagation of strangelets in 
the Earth's atmosphere, both of which are helpful to confirm their existence. Recently, 
Madsen has estimated the flux of strangelets in cosmic rays incident on the Earth [13J. 
As for the latter, unfortunately, the necessary interaction input physics is at best poorly 
known. 

Different phenomenological models of strangelet penetration in the atmosphere have 
been used by several authors to explain exotic cosmic ray events. Wilk et al. conjectured 
that although the initial mass of strangelets might be quite large, it decreases rapidly 
due to collisions with air molecules, until the mass reaches a critical value below which 
the strangelet disintegrates into neutrons [13] . Banerjee et al. provided a quite different 
scenario in which a strangelet picks up mass from atmospheric atoms |15j . Monreal 
novelly discussed the issue of strangelet accumulation in the atmosphere [T6] . 

In the present work, we will reinvestigate this issue within the framework of rotating 
liquid drop model, which is still a phenomenological model. We assume that SQM 
nuggets produced from any cosmological or astrophysical objects do reach the surface 
of the atmosphere, and we evaluate their behavior in the atmosphere. We find that 
strangelets with particular initial baryon number and particular initial Lorentz factor 
can reach mountain altitudes, even the sea level (~ 1000 g cm -2 ). 

In the following sections, we will first provide revised results of ground state 
strangelet calculated from liquid drop model. Then we will investigate the colliding cross 
section between a liquid strangelet and an air nucleus. Finally we will give numerical 
results about the propagation of strangelets in the atmosphere. It should be mentioned 
that, for the sake of simplicity, our calculation is limited to the ordinary (unpaired) 
strangelets, i.e. the possible color-superconductivity effect of strangelets has not been 
considered. 

2. Ground state properties of strangelets 

It was argued that the phenomenal bag model first used by Alcock & Farhi [17] may 
not suit for strangelets. Nevertheless, at the present level, bag model is still the most 
effective way to understand the properties of strangelets, such as ground state energy 
per baryon, charge-to-mass ratio, fissionability, etc. 

Within the framework of liquid drop model, He et al. [18J studied ground state 
properties of strangelets at finite temperature, by minimizing free energy of the system 
at fixed baryon number. The Coulomb energy was neglected there because the term 
contributes little to the system energy. In this section, we just go one step further to 
include Coulomb energy. Although Coulomb energy is negligible in computing E/A, it 
can greatly affect Z/A, and hence the fissionability parameters. 

We consider strangelet as gas of u,d,s quarks, their antiquarks, and gluons confined 
in an MIT bag model. The grand potential of the system Q = &i + BV, where B is 
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the bag constant, V is the volume, and the grand potential of species i is 

a = TT J q °° dk Pi (k) log(l ± exp(-(^/fc2 + m ? - Pi)/T)). (1) 

In the above equation, "±" denotes for Fermions/Bosons, \x is the chemical potential, 
p(k) denotes the density of states, which is given by 

where S (= AnR 2 for a sphere) is the surface area, C (= 87ri? for a sphere) is the 
curvature. The surface and curvature term for quarks are fg 9 \m/k) = — (l/87r)(l — 
(2/7r)arctan(/c/m)), f^ q \m/k) = (1/(12tt 2 ))(1 - (3k/2m)(n/2 - arctau(Jfe/m))), 
respectively; for gluons, = 0, = — 1/(Qtt 2 ). 
The free energy of the system is given by 

F = (fti + Nan) + E coul + BV, (3) 

i 

in which the term of Coulomb energy E cou i = (3/5)aZ 2 /R if electric charge is uniformly 
distributed in the sphere, where a is the fine structure constant and Z is the total 
electric charge. Since Coulomb energy is taken into account, the chemical potential of 
up quark and down quark (strange quark) are no longer identical. 

By minimizing F, the chemical potentials, charge to mass ratio and energy per 
baryon at any given baryon number and temperature can be calculated. Taking strange 
quark mass m s = 150 MeV and bag constant B = (145 MeV) 4 , we get the following 
fitting values for the parameters. At zero temperature, 

M str /A = (314.6(4)A- a532(4) + 875.9(1)) MeV, (4) 

therefore, the minimum baryon number for stability (Mj A < 930 MeV) is A crit = 27. 
The surface energy is 

E s /A 2/3 = (69.0(2)A-°' 466(3) + 77.9(1)) MeV, (5) 

and the rescaled radius is 

r = R/A 1/3 = (0.124(l)A- a445(3) + 0.941(1)) fm, (6) 

The numbers in parenthesis following each value indicate just the fitting uncertainties 
of the value in the last digit. The charge-to- mass ratio with respect to the equivalent 
baryon number is shown in Fig.fCQ). 

As for a strangelet at excitation states to de-excite, 7— ray emission, hadron 
emission and fission into small parts should be under consideration. 7— ray emission 
and meson emission do not change the baryon number a little. According to CEFT 
model [19], baryon evaporation is suppressed in term of meson evaporation due to 
much smaller probability to simultaneously form two quark- ant iquark pairs than one 
pair. Banerjee et al. [20] and Sumiyoshi et al. [21] has calculated meson and baryon 
evaporation rate of QGP (p q = 0), respectively. As for a strangelet, in which the quarks 
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Figure 1. Charge to mass ratio of strangelets. 



have non-zero chemical potentials, the numerical results are as follows. The energy loss 
rate caused by baryon evaporation is 
dE, 



dt 



\baryon = -6.23 x 10 19 A 2/3 T 2 exp(-999.9/T) MeV sec" 1 , 



(7) 



(8) 



while the energy loss rate caused by meson evaporation is 

fU,^-, 12Xl W e Xp (-3e, 1/r )MeV S ee-V 

Therefore, in the de-excitation process the only effective way to change the baryon 
number of the strangelet is the fission. Note that a strangelet at excited states will 
rapidly release its energy in about 10 -18 ~ 10~ 15 s, which is negligible compared with 
colliding time intervals. 



3. Fission of strangelets by colliding with air nuclei 

Now we investigate the stability of strangelets using the rotating liquid drop model. 
In the case of non-rotating systems, the relation between the nature of the stationary 
points and the stability of a system is simple, a maximum in one or more degrees of 
freedom indicates instability. However, the case for rotating systems is more subtle. 

Consider a configuration of a rotating incompressible uniformly charged fluid 
endowed with a surface tension. The effective potential energy is given by 

E = E s + E c + E R , (9) 

where Eg is the surface energy, Ec the electrostatic energy, and Er the rotational 
energy. We neglect the curvature energy because it contributes little to the issue we 
consider. 
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We may write the deformation energy, measured with respect to the energy of the 
sphere, in the following dimensionless form, familiar in the literature of nuclear fission, 

E - E<® 

£ = 755— = (B s - 1) + 2x(Bc - 1) + y(B R - 1). (10) 

E s 

Here B s = E s /E s °\ B c = E c /E§\ B R = E R /E$\ and the two dimensionless 
parameters x and y specify the ratios of electrostatic and rotational energies of the 
sphere to the surface energy of the sphere, which are defined 4 0) /24 0) and 

y = Eft /Eg . According to our calculations, we found that the fissionability parameter 
x of strangelets varies from 0.001 to 0.030, which is much smaller than normal nuclei. 
Therefore, we take x = for simplicity in the following calculations. 

If there is no rotation (y — 0), the ground state is a sphere and the saddle shape 
is the configuration of two tangent spheres. With increased rotation, the ground state 
sphere is flatten into an axially symmetric (Hiskes) shapes, and the saddle varies to 
the so-called Pik-Pichak shapes. If y is even larger, the ground state will convert to 
a triaxial (Beringer-Knox) shape, which is quite close in appearance to the Pik-Pichak 
saddle. In fact, there exist a critical value y cr a above which the fission barrier vanishes, 
which implies a maximum rotation for stability. 

Cohen et al. [22] has calculated the ground state and fission barrier energy measured 
with respect to the energy of the sphere, 

Cground = y(-0.056 + 0.04% - 1.35% 2 + 0.946y 3 ), (11) 
Carrier = 0.280 - 0.77% + 0.622?/ 2 - 0.105?/ 3 . (12) 

For x = 0, y crit = 0.79. 

Now we consider a cosmic ray strangelet incident in the Earth's atmosphere. The 
center-of-mass energy E cm = (M 2 r + M 2 air + 2^M str M air ) 1 ^ - M str - M air , where 
M a i r = 14M (M : approximately the proton mass) is the mass of the air nucleus, 
M str and 7 is the mass and the Lorentz factor of the strangelet, respectively. There 
exists a nonzero colliding cross section as long as the strangelet has enough kinetic 
energy to overcome the Coulomb barrier E v between the strangelet and the air nucleus, 
i.e. if E cm > E v , they will have a chance to "fusion" into a compound strangelet. In 
other words, the strangelet "absorbs" the air nucleus. 

After fusion of the two, the excitation energy E e = E cm — E v — E^' — E groun d. 
According to rotating liquid drop model, if the projectile strangelet is energetic enough, 
E e will be higher than the fission barrier Eb, the newly formed compound strangelet 
will fission into two smaller strangelets which have nearly equal baryon numbers. 

However, with the increasing kinetic energy of projectile strangelet, the effect of 
rotation should no longer be neglected, since if the fission barrier of rotating compound 
system approaches zero, no compound strangelet will form. It is reasonable to suppose 
the interaction time scale in this case is much shorter than the former. 

The geometric cross section for contact is a geo = tt^qA]^ + 1.12A^{ 3 ) 2 , which is 
used by some earlier studies for crude calculations. However, the cross section is not 
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Figure 2. Critical value of strangelet velocity. The upper curve represents the fission 
threshold of a non-rotating strangelet, and the lower curve represents the collision 
threshold. 

always o geo for different baryon numbers and different velocities. The cross section for 
close collision can be put as 

&col = 0geo(l Z s t r Z a i r e / (R s tr + Rair)E cm ) , (13) 

and the cross section for fusion is 

a cn = 2vy cnt I Ef\M 2 str + M 2 air + 2 1 M str M air )/(( 1 2 - l)M 2 tr M 2 air ). (14) 

Fig. (El) gives the critical value of velocity for fusion and fission as a function of 
baryon number. The relation between the impact parameter b in a collision between the 
two masses and the velocity of the strangelet is shown in Fig. (J3j). 

4. Propagation of strangelets in the atmosphere 

We consider strangelets with zero zenith angle of trajectory in this work, and the 
influence of gravity force is neglected because of its small contribution to the issue. 
Therefore, strangelet-air collision and ionization effect are what we mainly concern. 

In our model, there exists a critical velocity of the strangelet in the strangelet-air 
collisions below which the air nucleus will be fused with the strangelet, and above which 
some mass will be stripped from the strangelet. The model is based on the analogical 
result in nuclear collisions, i.e. the linear momentum transfer between the strangelet 
and the air nucleus reaches a maximum around some critical energy. 

When fusion happens, the velocity of the strangelet drops to 

7 ' = {lM str + M air )/(M 2 tr + M 2 ir + 2 1 M stT M air ) 1 ' 2 (15) 
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Figure 3. A collision diagram of square of the impact parameter b (times it) versus 
the velocity of a strangelet with A str = 100, for the bombardment of an air nucleus by 
the strangelet. 

after each collision, and particularly if E e > E{,, the strangelet will fission into two 
smaller strangelets with equal baryon number 0.5(A str + A air ) and probably equal 
longitudinal velocity. 

When the strangelet is more energetic, no compound strangelet will form. If 
the experimental law in nuclear collision is adaptable in this case, the velocity of the 
strangelet is assumed to drop to 

7 ' = 7 - {icri -Tcrt)> ( 16 ) 

where ^ cri is the critical gamma factor for fusion, and ^' cri can be deduced from Eq lfToT) 
given 7 = 7 cri . The value of 7 cr j can be found by solving the equation a cri = 0.5 <j co i, i.e. 
the watershed of fusion-dominated and stripping-dominated collisions. It should also be 
mentioned that, we assume new strangelet will have a baryon number of (A str — A air ) 
after each collision in numerical calculations. Indeed, at the present level, the mass 
and energy spectrum of decay products in this range are quite uncertain. Although our 
supposition is somewhat crude, it nevertheless tells some important information. 

In addition to the effect of colliding with air nuclei, the issue of the energy loss 
of the strangelet through ionization of surrounding media can not be avoided, which is 
described by the Bethe-Bloch stopping power formula [23] . 

dE/dx = -0.153/r 2 Z s 2 tr (ln(7 2 - 1) - /3 2 + 9.39) MeV/(g cm" 2 ) (17) 

If v < voZglr , where vq = 2.2 x 10 8 cm/s is the speed of electron in the first Bohr orbit, 
the effective charge Z str — > (v/v )Zj r due to the effect of electron capture [24J. 
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The following figures show our numerical results. In Fig. (j5J) , we present the mass 
evolvement of strangelets with 7 = 10 3 as a function of atmospheric depth. It is obvious 
that the larger the 70 and Aq, the more possible the strangelet reach the sea level. In 
Fig-©, we show the distribution of particular final (x = 1036 g cm~ 2 ) baryon number of 
strangelets as a function of initial baryon number and gamma factor, the upper serried 
oblique lines correspond to products of fission (lower energy), and the lower transverse 
lines correspond to products of stripping (higher energy). We find that if the initial 
strangelets have A > 3000 or 70 > 140, they will have a chance to be detected by 
ground-based experiments. In Fig. ([6]), we give the final Lorentz factor as a function of 
initial baryon number. 
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Figure 4. Mass evolvement of strangelets with A = 1000, 3000, 5000, 10000 and 
70 = 10 3 as a function of atmospheric depth. 



5. Conclusions 

From above discussion, it is reasonable to make a conclusion that strangelets with large 
mass and energy have the chance to penetrate the atmosphere to reach the sea level. 
Our model gives a lower limit of initial baryon number, that is, A > 3000 or 70 > 140. 
Relevant flux for ordinary strangelets is unclear yet, but as predicted by Madsen [IB], it is 
about 1 ~ 10 per sqm year sterad, so there is a possibility for ground-based experiments 
to detect them. Madsen [13] predicts a flux of strangelets with a velocity spectrum (an 
event with 7 = 140 could be very unlikely there) in a model where strangelets originate 
only by merging of binary strange stars. However, the 7 factor of strangelets produced 
in other ways (e.g., ejected and accelerated in pulsar's magnetospheres [HI EH], [26]) may 
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Figure 5. Final mass distribution of strangelets 
initial baryon number and lorentz factor. 



1036 g cm 2 ) as a function of 



100 




3000 4000 5000 6000 7000 8000 
initial baryon number^ 



9000 10000 



Figure 6. Final lorentz factor (x — 1036 g cm 2 ) 
number with 70 = 200, 400, 600, 800, 1000. 



as a function of initial baryon 



be greater than 140, and the possibility of an event with higher 7 could then not be 
ruled out yet. 
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